!****************************************************
! Default cosmological parameters, used by object.F90
!****************************************************

module mParameters

        use mPrecision

#ifndef PARAMETER_SET
        real(D), parameter :: DEFAULT_OMEGA_M0 = 0.310_D, &
                              DEFAULT_OMEGA_V0 = 0.690_D, &
                              DEFAULT_W_0 = -1.0_D, &
                              DEFAULT_W_A = 0.0_D, &
                              DEFAULT_OMEGA_B0 = 0.0528_D, &
                              DEFAULT_H = 0.660_D, &
                              DEFAULT_SIGMA_8 = 0.882_D, &
                              DEFAULT_B_1 = 0.998_D, & ! lowest-order galaxy bias
                              DEFAULT_B_2 = 0.0_D, &   ! next-order galaxy bias
                              DEFAULT_MASS_NU = 0.0_D, &
                              DEFAULT_SPECTRAL_INDEX = 0.979_D, &
                              DEFAULT_RUN = 0.0_D, &
                              DEFAULT_A_S = 0.82 ! CMB normalization parameter

#elif PARAMETER_SET == 2
        real(D), parameter :: DEFAULT_OMEGA_M0 = 0.27_D, &
                              DEFAULT_OMEGA_V0 = 0.73_D, &
                              DEFAULT_W_0 = -1.0_D, &
                              DEFAULT_W_A = 0.0_D, &
                              DEFAULT_OMEGA_B0 = 0.046_D, &
                              DEFAULT_H = 0.72_D, &
                              DEFAULT_SIGMA_8 = 0.9_D, &
                              DEFAULT_B_1 = 1.0_D, & ! lowest-order galaxy bias
                              DEFAULT_B_2 = 0.0_D, & ! next-order galaxy bias
                              DEFAULT_MASS_NU = 0.0_D, &
                              DEFAULT_SPECTRAL_INDEX = 1.0_D, &
                              DEFAULT_RUN = 0.0_D, &
                              DEFAULT_A_S = 0.87 ! CMB normalization parameter

#elif PARAMETER_SET == 1
        real(D), parameter :: DEFAULT_OMEGA_M0 = 0.3_D, &
                              DEFAULT_OMEGA_V0 = 0.7_D, &
                              DEFAULT_W_0 = -1.0_D, &
                              DEFAULT_W_A = 0.0_D, &
                              DEFAULT_OMEGA_B0 = 0.05_D, &
                              DEFAULT_H = 0.7_D, &
                              DEFAULT_SIGMA_8 = 0.9_D, &
                              DEFAULT_B_1 = 1.0_D, & ! lowest-order galaxy bias
                              DEFAULT_B_2 = 0.0_D, & ! next-order galaxy bias
                              DEFAULT_MASS_NU = 0.0_D, &
                              DEFAULT_SPECTRAL_INDEX = 1.0_D, &
                              DEFAULT_RUN = 0.0_D, &
                              DEFAULT_A_S = 0.87 ! CMB normalization parameter

#endif ! PARAMETER_SET == 1

        public

end module mParameters
